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Abstract. We present recent results on light mesons based on Dalitz plot analyses of charm decays 
from Fermilab experiment E791. Scalar mesons are found to have large contributions to the decays 
studied, D + — > K~K + K + and D + ,D+ — > K~K + K + . From the Kkk final state, we find good evidence 
for the existence of the light and broad K meson and we measure its mass and width. We also discuss 
recently published results on the 3ji final states, especially the measurement of the /o parameters 
and the evidence for the o meson from D + — » an + . These results demonstrate the importance of 
charm decays as a new environment for the study of light meson physics. 



INTRODUCTION 

The decays of charm mesons are currently a new source of information for the study of 
light meson spectroscopy, with the advantages of having well defined initial state (the 
D meson, a state with defined mass). This new information is complementary to that 
from scattering experiments and can be particularly relevant to the understanding of the 
scalar sector. 

Here we present preliminary results for the Dalitz-plot analysis of the Cabibbo- 
favored decay D + —> K~k + k + using data from Fermilab E791 experiment. We also 
present an overview of our results for D+ — > 7T.~7r. + 7T. + [1] and D + — » 7i~7i + 7i + [2] 
Dalitz-plot analyses. The E791 data was collected in 1991/92 from 500 GeV/c 7tT- 
nucleon interactions. For details see [3]. 

For the D + —> K~n + n + analysis, when we include all known Kn resonant channels 
plus a non-resonant (NR) contribution, we find that the NR decay is dominant. This is 
unusual in D decays. Moreover, the fit model has important discrepancies with respect 
to the data. By including an extra scalar resonant state, with unconstrained mass and 
width, we obtain a fit which is substantially superior to that without this state. The 
values for its mass and width are found to be 797 ± 19 ± 42 MeV/c 2 and 410 ± 43 ± 85 
MeV/c 2 respectively. We refer to this state as the K. The existence of such a state has 
been greatly discussed in the literature in recent years [4]-[14]. We also obtain new 
measurements for the mass and the width of the ATq(1430) resonance: 1459 ±7 ±6 
MeV/c 2 and 175 ± 12 ± 12 MeV/c 2 respectively. 

From our analysis of D+ — > 7T.~7r. + 7T + decays, we find that the dominant decay fraction 



1 For the E791 Collaboration. 



comes from /o(980)7t + . We obtain new measurements for the /o(980) and /o(1370) 
masses and widths. From D + — > 7t~7t + 7t + decays, we find that a model with only known 
nn resonances plus a NR channel is not able to describe the data adequately. We find 
strong evidence for the presence of a light and broad scalar resonance, the o(500), the 
on + channel being responsible for half of the decay rate. We measure the mass and 
the width of this scalar meson to be 478+^ ± 17 MeV/c 2 and ±21 MeV/c 2 , 

respectively. 



THE D + -> K~% + % + DALITZ-PLOT ANALYSIS 

From the original 2 x 10 9 events collected by E791, and after reconstruction and selec- 
tion criteria, we obtained the D + — > K~k + k + sample shown in Figure 1(a). The filled 
area represents the level of background; besides the combinatorial, the other main source 
of background comes from the reflection of the decay Df — > K K + n + (through K*K + 
and The crosshatched region contains the events selected for the Dalitz-plot ana- 
lysis. There are 15090 events in this sample, of which 6% are background. 

Figure 1(b) shows the Dalitz-plot for these events. The two axes are the squared 
invariant-mass combinations for K%, and the plot is symmetrized with respect to the 
two identical pions. The plot presents a rich structure, where we can observe the clear 
bands from K* (890)7t + , and an accumulation of events at the upper edge of the diagonal, 
due to heavier resonances. To study the resonant substructure, we perform an unbinned 
maximun-likelihood fit to the data, with probability distribution funtions (PDF's) for 
both signal and background sources. In particular, for each candidate event, the signal 
PDF is written as the square of the total physical amplitude jz (defined below) and it 
is weighted for the acceptance across the Dalitz plot (obtained by Monte Carlo (MC)) 
and by the level of signal to background for each event, as given by the line shape of 
Figure 1(a). The background PDF's (levels and shapes) are fixed for the Dalitz-plot fit, 
according to MC and data studies. 

We begin describing our first approach to fit the data, which represents the conven- 
tional Dalitz-plot analysis including the known Kit resonant amplitudes (A n , n > 1), plus 
a constant non-resonant contribution. The signal amplitude is constructed as a coherent 
sum of the various sub-channels: 

N 
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Each resonant amplitude is written as 

Sl n = BW n F^ F<! J) rti J) . (2) 
where BW n is the relativistic Breit-Wigner propagator 

BW n = — ; 1 , — - (3) 



m. 



q — m — imoT(m) 



FIGURE 1. (a) The K % + % + invariant mass spectrum. The filled area is background; (b) Dalitz plot 
corresponding to the events in the dashed area of (a). 



with mass dependent width, 

The quantities Fp and Fr are the Blatt-Weisskopf damping factors respectively for the 
D and the Kit resonances, they depend on the radii of the decaying meson and are set to 
r D = 3.0 GeV 1 and r R = 1.5 GeV 1 [18]; p* is the pion momentum in the resonance 

rest frame at mass m\2 {p* Q = p*(mo)). M,^ describes the angular distribution due to the 
spin / of the resonance. See details in [2]. Finally each amplitude is Bose symmetrized 
^ = ^[(12)3]+^„[(13)2]. 

Using this model (Model A), we find contributions from the following channels: 
the non-resonant, responsible for more than 90% of the decay rate, followed by 
£ *( 1430)71+, iT (892)7t+, K*( 1680)tt+ and £|(1430)ji + . The decay fractions and rela- 
tive phases are shown in Table 1 . These values are in accordance with previous results 
from E691 [15] and E687 [16]. We thus confirm a high non-resonant contribution 
according to this model, which is totally unusual in D decays. Besides, there is an 
important destructive interference pattern, since all fractions add up to 140 %. 

To evaluate the fit quality of Model A, we compute a two-dimensional % 2 in the 
Dalitz plot, from the difference in densities for the model (from a fast-MC algorithm) 
and the data. We obtain % 2 /v = 2.7 (v being the number of degrees of freedom), 
with a corresponding confidence level (CL) of 10~ H . In Figure 2(a) we show the Kn 
low and high squared-mass projections for data (error bars) and model (solid line). 
The discrepancies are evident at the very low-mass region for m 2 (Kn)i ow and near 



TABLE 1. Results without K (Model A) and with K (Model B). Preliminary. 



Decay 
Mode 


Model A: No K 
Fraction (%) Phase 


Model B: With K 
Fraction (%) Phase 


NR 

K7C+ 

£*(892)tc+ 
£ *( 1430)71+ 
£ 2 * (1430)71+ 
£*(1680)tc+ 


90.9 ±2.6 0° (fixed) 

13.8±0.5 (54±2)° 
30.6 ±1.6 (54 ±2)° 
0.4±0.1 (33±8)° 
3.2±0.3 (66±3)° 


13.0±5.8±2.6 (349±14±8)° 
47.8±12.1±3.7 (187±8±17)° 
12.3 ±1.0 ±0.9 0° (fixed) 
12.5±1.4±0.4 (48±7±10)° 
0.5±0.1±0.2 (306±8±6)° 
2.5 ±0.7 ±0.2 (28 ±13 ±15)° 



2.5 (GeV/c 2 ) 2 for m 2 (Kn) high . These regions of disagreement are the same observed 
previously by E687 [16]. We thus conclude that a model with the known Kn resonances, 
plus a non-resonant amplitude, is not able to describe the D + — > K ~7T. + 7T. + Dalitz plot 
satisfactorily. 

A similar pattern - bad fit quality with large NR fraction - is found in the analysis of 
the decay D + — > 7T.~7r. + 7r. + when allowing only the established 7UT. resonances [2] . There 
we find that the inclusion of an extra scalar resonance improves the fit substantially, 
giving strong evidence for the c(500). See the section on D + — > 7T.~7r. + 7r. + below. Thus, 
we are lead to try an extra scalar resonance in our fit model here. The possible existence 
of a light and broad Kit scalar state has been suggested by many authors [4, 5, 6, 7, 8], 
some of them believing it would be a member of a light scalar nonet [9, 10, 1 1]; however, 
its existence has been the subject of some controversy also[12, 13, 14]. 

A second fit model, Model B, is constructed by the inclusion of an extra scalar 
state, with unconstrained mass and width. For consistency, the mass and width of the 
other scalar state, the £q(1430), are also free parameters of the fit. We adopt a better 
description for these scalar states by introducing gaussian-type form-factors [12] to take 
into account the finite size of the decaying mesons. Two extra floating parameters are 
the meson radii r D and r R introduced above. 

Using this model, we obtain the values of 797 ± 19 ± 42 MeV/c 2 for the mass and 
410 ± 43 ± 85 MeV/c 2 for the width of the new scalar state (first error statistical, second 
error systematic), referred to here as the K. The values of mass and width obtained for the 
*o (1430) are respectively 1459 ± 7 ± 6 MeV/c 2 and 175 ± 12 ± 12 MeV/c 2 , appearing 
heavier and narrower than presented by the PDG [17]. The decay fractions and relative 
phases for Model B, with systematic errors, are given in Table 1 . Compared to the results 
of Model A (without k), the non-resonant mode drops from over 90% to 13%. The K7t + 
state is now the dominant channel with about 50%. The meson radii ro and rR are found 
to be respectively 5.0±0.5 GeV 1 and 1.6± 1.3 GeV \ in complete agreement with 
previous estimates [18, 19]. 

Moreover, the fit quality of Model B is substantially superior to that of Model A. The 
% 2 /v is now 0.73 with a CL of 95%. The very good agreement between the model and 
the data can be seen in the projections of Figure 2(b). 

A number of studies were done to check these results. For example, we replaced the 
complex K Breit-Wigner by a real Breit-Wigner, with no phase variation. In this case, 
we got similar mass and width for this extra state, but with unphysical fractions for 
this state and the NR, and a worse fit quality. We also replaced the K by a hypothetical 




vector state, with unconstrained mass and width, but it appears with a small fraction, 
the fit quality being comparable to the model without it, and all fractions and phases 
remaining unchanged. A tensor model was also tried without convergence, the width 
being driven to large, negative values. Other models with the K were also tried. For 
example, modifications to the scalar Breit-Wigner amplitude and to the form-factors 
were introduced. A number of studies for the parameterization of the NR amplitude 
were tried [22], with and without the K. No model without the K was able to describe our 
data satisfactorily. All variations of models with K gave similar results for the K mass 
and width (within errors) although the fractions for K7t and NR showed correlations. 

Thus, from the results above, we find strong evidence that a light and broad scalar Kit 
resonance gives an important contribution to the D + — > K~% + % + decay. 



THE D+ -> 71 71+71+ RESULTS 

In Figure 3 we show the invariant mass distribution for the sample collected 

by E791 after reconstruction and selection criteria [1, 2]. Besides combinatorial back- 
ground, reflections from the decays D + — > K~n + ii + , D + — > K n + (plus one extra 
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FIGURE 4. (a) The Df -> 7i~7i + 7i + Dalitz plot and (b) the D + -> 7i~7i + 7i + Dalitz plot. Since there are 
two identical pions, the plots are symmetrized. 



track) and Df — > r| / 7T + , r\ f — » p°(770)y are all taken into account. The hatched re- 
gions in Figure 3 show the samples used for the Dalitz-plot analyses. There are 1686 
and 937 candidate events for Z) + and D+respectively, with a signal to background ratio 
of about 2:1. The Dalitz plots for these events are shown in Figure 4, the axes corre- 
sponding to the two invariant-masses squared. For the Dalitz-plot fits of both 
Z) + —¥ 7l~7l + 7i: + and Df —* 7C~7l + 7l + decays, we use essentially the same formalism as 
for the D + — > K~n + n + decays. See details in [1, 2]. 



TABLE 2. Dalitz fit results for D+ -> 7C-7C + 7C + . 



Decay Mode 


Fraction (%) Phase 


/o(980)7t+ 

NR 
p°(770)7t+ 
£(1270)71+ 
/o(1370)ji+ 
p°(1450)7t+ 


56.5 ±4.3 ±4.7 0° (fixed) 
0.5±1.4±1.7 (181±94±51)° 
5.8±2.3±3.7 (109±24±5)° 
19.7±3.3±0.6 (133±13±28)° 
32.4±7.7±1.9 (198±19±27)° 
4.4±2.1±0.2 (162±26±17)° 



For the D+ — > % events in Figure 4(a), the signal amplitude includes the 

following channels: p°(770)7r+, / (980)ji+ / 2 (1270)ji+ / ( 1370)71+, p°(1450)tt + and 
the non-resonant, assumed constant across the Dalitz plot. 

For the /o(980)7t + amplitude, instead of a simple Breit-Wigner of Eq. 3 2 , we use a 
coupled-channel Breit-Wigner function [23], 

r n = Sn\jmlJA-ml, ^k = y (^JmlJ4-m 2 K+ + \JmlJA-m 2 K ^j . (6) 

The Df — > k~k + k + Dalitz plot is fit to obtain not only the decay fractions and phases of 
the possible sub-channels, but also the parameters of the /o(980) state, g n , gK, and mo, 
as well as the mass and width of the /o(1370). The other resonance masses and widths 
are taken from the PDG[17]. The resulting fractions and phases are given in Table 2. 

The measured / (980) parameters are m = 977 ± 3 ± 2 MeV/c 2 , g n = 0.09 ± 0.01 
± 0.01 and gK = 0.02 ± 0.04 ± 0.03. Our value for g n is in very good agreement with 
OPAL and MARKII results [24], but WA76 [23] found a much larger value, g n = 0.28 
± 0.04. Our value of gK indicates a small coupling of /o(980) to KK. The values of the 
/o(980) mass and of g n , as well as the magnitudes and phases of the resonant amplitudes, 
are relatively insensitive to the value of gK- Both OPAL and MARKII results are also 
insensitive to the value of gK. WA76, on the contrary, measured gK = 0.56 ± 0.18. 

By fitting the Dalitz plot using for the /o(980) a simple Breit-Wigner function, we 
find m = 975 ± 3 MeV/c 2 and T = 44 ± 2 ± 2 MeV/c 2 , and the results for fractions and 
phases are indistinguishable. 

The confidence level of the fit for Df — > 7T.~7r. + 7r. + is 35% [1]. In Figure 5 we show the 
7T.~7T. + mass-squared projections for data (points) and model (solid lines, from fast-MC). 

As we can see by the results of Table 2, approximately half of the Df — > 7T.~7r. + 7r. + rate 
is via /o(980)7T. + . If the spectator amplitude is dominant in this decay, this would support 
the interpretation of the /o(980) as an ss state. On the other hand, the large contribution 
from /o(1370)7T. + indicates the presence of either W-annihilation amplitudes or strong 
rescattering in the final state. In fact, the /o(1370)7l + is not observed in the D+ — > 



2 For both D,D S — ► 371 analyses, the relativistic Breit-Wigner for each resonant amplitude is defined with 
a factor (-1) with respect to Eq. 3. 
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FIGURE 5. 5i2 and 513 (m^ K ) projections for D+ — ► % 7t + 7t + data (dots) and our best fit (solid). The 
hashed area corresponds to background. 



K + K % + final state[25], pointing to the /o(1370) being a non-55 state, as suggested 
by the naive quark model[17]. There is no evidence in the D+ decay for a low-mass 
broad scalar particle as seen in the D + — >■ n~n + 'K, + decay, discussed below. 



THE D+ -> 7171+71+ DECAY 

In a first approach, we try to fit the D + — > k~k + k + Dalitz plot of Figure 4(b) with 
the same amplitudes used for the Df — > k~k + k + analysis. Using this model, the non- 
resonant, the p°(1450)7T. + , and the p°(770)7l+ amplitudes are found to dominate, as 
shown in Table 3, and in agreement with previous reported analyses [26, 27]. However, 
this model does not describe the data satisfactorily, especially at low k~k + mass squared, 
as can be seen from Fig. 6(a). The % 2 /v obtained from the binned Dalitz plot for this 
model is 1.6, with a CL less than 10~ 5 . 

To investigate the possibility that another 7171+ resonance contributes to the D + — > 
7r.~7r. + 7T. + decay, we add an extra scalar resonance amplitude to the signal PDF, with 
mass and width as floating parameters in the fit. 

We find that this model improves our fit substantially. The mass and the width of 
the extra scalar state are found to be 478+23 + 17 MeV/c 2 and 324+ ±21 MeV/c 2 , 
respectively. Refering to this state as the o, we obtain that the ott + chanell produces 
the largest decay fraction, as shown in Table 3; the non-resonant amplitude, which is 
dominant in the model without ott + , drops substantially. This model describes the data 
much better, as can be seen by the 7T.Tr. mass squared projection in Fig. 6(b). The % 2 /v is 
now 0.9, with a corresponding confidence level of 91%. 

The existence of a light tutc state, or the o, has been the subject of a long-standing 
controversy [28, 29]. Various experiments have presented inconsistent evidence for this 
state [30], yielding conflicting results [17, 29]. 




TABLE 3. Dalitz fit results for D + — > % 7C+7C+ . First errors are statistical, second 
systematics (only for fit with OTC+ mode). 



Decay 
Mode 


Fit without C7C+ 
Fraction (%) Phase 


Fit with <J7C+ 
Fraction (%) Phase 


crct+ 
P (770)tc+ 

NR 
/o(980)7t+ 
f 2 (1270)71+ 
/o (1370)tc+ 
P°(1450)jc+ 


20.8 ±2.4 0° (fixed) 
38.6±9.7 (150±12)° 
7.4±1.4 (152±16)° 
6.3 ±1.9 (103 ±16)° 
10.7±3.1 (143±10)° 
22.6 ±3.7 (46 ±15)° 


46.3±9.0±2.1 (206±8±5)° 
33.6 ±3.2 ±2.2 0° (fixed) 
7.8±6.0±2.7 (57±20±6)° 
6.2±1.3±0.4 (165±11±3)° 
19.4 ±2.5 ±0.4 (57 ±8 ±3)° 
2.3±1.5±0.8 (105±18±1)° 
0.7±0.7±0.3 (319±39±11)° 




FIGURE 6. s\2 and su (m^ K ) projections for D + — > 7C~7C+7C+ data (dots) and our best fit (solid) for 
models (a) without and (b) with OTC+ amplitude. The dashed distribution corresponds to the expected 
background level. 



To test the model above, we replace the scalar amplitude by vector and tensor states, 
and also by a real Breit-Wigner, with no phase variation (as also done in the D + — > 
K~n + n + analysis). All these alternative models fail to describe the data as well as the 
scalar (regular) Breit-Wigner amplitude. See detailed discussion in [2]. 



CONCLUSION 

From the data of the Fermilab E791 experiment, we studied the Dalitz plots of the decays 
D + — > K~k + k + , D+ — >• 71^71+71+ and D + — > k~k + ti + . In these three final states, the 
scalar intermediate resonances were found to give the main contribution to the decay 
rates. We obtained strong evidence for the existence of the a and k scalar mesons, 
measuring their masses and widths. We also obtained new measurements for masses 
and widths of the other scalars studied, / (980), /o(1430) and ^(1430). 

The results presented here show the potential of D meson decays for the study of light 
meson espectroscopy, in particular in the scalar sector. 
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